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I. INTRODUCTION 

Onsager's intuition that purely repulsive rods undergo 
an entropy-driven transition from an isotropic (/) to an 
orientationally ordered nematic [N] phase constitutes 
one of the major milestones in our understanding of liquid 
crystals [T] . The key ingredient of this phenomenon relies 
on considering markedly non-spherical particles, which 
can be modeled as cylindrically symmetric "rods" and 
"plates". In the early 1970s Preiser pointed out that a 
richer phase behavior is expected, if the assumption of 
cylindrical symmetry is released [51. Besides the usual 
prolate (Af+) and oblate {N-) uniaxial nematic phases, 
normally developed by uniaxial rods and plates, respec- 
tively, a novel nematic phase with an increased orienta- 
tional order can appear in the phase diagram. Such a 
liquid crystal phase is characterized by alignment along 
three directors and, consequently, by the presence of two 
distinct optical axes, hence the name biaxial nematic 
[Nb) [3]. Further studies suggested that Nb stability 
could be interpreted as a balanced competition between 
rodlike (favoring and platelike (favoring N^) behav- 
ior [4-6J. 

In more than 40 years since its first theoretical predic- 
tion, extensive theoretical [MTTj and simulation 
work has been devoted to identify the conditions under 
which a stable Nb phase could be observed. The practi- 
cal limitations in this sense are testified by the fact that, 
apart from the micellar system studied by Yu and Saupe 
[25] . no such state has been observed for more than 30 
years. A renewed interest towards the topic has grown 
due to the first experimental realization of thermotropic 
Nb liquid crystals in systems of bent-core molecules a few 
years ago [26l[22]. In lyotropics, a remarkably stable Nb 
phase was recently discovered in a colloidal suspension of 
mineral boardlike particles 

Boardlike particles, that is, particles with the symme- 



try of a brick, represent the simplest model in which an 
Nb phase has been predicted [3]. However, the emer- 
gence of smectic layering is expected to prevent the real- 
ization of this phase, unless the constituent particles are 
designed with a precision far beyond present-day ability 
[H [2H] • A higher Nb stability can be achieved by consid- 
ering size-polydisperse systems of boardlike particles, as 
demonstrated by a recent experiment 28]. In fact, poly- 
dispersity seems to enhance Nb stability through two 
distinct phenomena: (i) a reduced smectic stability |T2] 
and (ii) an N^-N- competition, which manifests itself ex- 
clusively in systems of slightly elongated (rodlike) boards 
|17) . The first phenomenon does not come as a surprise 
I12[ , since it is well known that polydispersity renders the 
establishment of long-distance positional ordering unfa- 
vorable [291 - I3T] . On the contrary, the reason behind the 
second phenomenon appears to be more obscure. 

In this paper we investigate the effect of a non- 
adsorbing depletant on the biaxial nematic stability of 
(monodisperse) boardlike particles. Our understand- 
ing of depletion dates back to the pioneering work by 
Asakura and Osawa [51] and Vrij [33], who showed that 
the addition of small co-solutes (e.g. polymers, sur- 
factants, micelles) to a colloidal suspension gives rise 
to an effective attraction between colloidal particles. 
Since then, the concepts related to depletion have been 
widely applied to various scientific fields [34]: in biology 
by interpreting phenomena like macromolecular crowd- 
ing ^35j and protein crystallization [36) : in nanotechnol- 
ogy through e.g. the development of self-assembly pro- 
cesses as key- lock structures [37l|38|; in condensed matter 
physics, furnishing answers to fundamental problems like 
the condition for gas-liquid phase separation [39] , the ki- 
netics of crystallization [dO] |4T] and the nature of glassy 
states [42] . More recently, the liquid crystal phase behav- 
ior of non-spherical colloids, typically rods [i3Mi8] and 
plates [I51 - f52] . in presence of a depletant has also been 



addressed. As a general feature, the addition of a deple- 
tant reduces the stabihty of hquid-crystal phases, leading 
to a direct isotropic-crystal transition at high enough de- 
pletant mole fraction. Moreover, when the size of the de- 
pletant particles is big enough, one or more critical points 
appear in the phase diagram, indicating a liquid-gas sep- 
aration between phases with same spatial symmetries. 

In contrast to the aforementioned work on rods and 
plates, we focus here on the low depletant density limit, 
where the stability of the nematic liquid crystal phases 
developed by the pure system of boardlike particles is 
preserved. In the same spirit as the Asakura-Oosawa- 
Vrij model for spheres [32 [33] , we consider the limit of 
low depletant density and neglect depletant-depletant in- 
teractions. For the sake of convenience, we model the 
depletants as cubic particles excluded from the surface 
of the cuboids via a hard-core interaction. A mean-field 
theory at second virial order [TJ [53] with restricted orien- 
tations (Zwanzig model) [51] constitutes our theoretical 
framework. The degrees of freedom of the depletant in 
the partition function can be systematically integrated 
out, giving rise to an effective potential between board- 
like particles [S51 [S3] , where only two-body interactions 
are considered. The assumption of ideal depletant allows 
to determine an explicit expression for such a pairwise 
depletion potential. We show that, by varying the deple- 
tant density, the system develops an N^-N^ competition 
remarkably similar to that predicted for a polydisperse 
system of boardlike particles in absence of depletant [T7] . 
If in Ref. [17j the origin of this competition is not ev- 
ident, here it appears to be due to a balance between 
the hard-core repulsion between boardlike particles, fa- 
voring ordering, and the depletion attraction, favor- 
ing N_ ordering. As a consequence of this effect, the 
biaxial nematic phase appears to be stable over a wide 
range of depletant density. We therefore suggest that the 
concentration of a non-adsorbing depletant furnishes in 
practical situations the simplest, though effective, way 
to control the liquid-crystal phase behavior of boardlike 
particles and to select states of high biaxial-nematic sta- 
bility. 

The paper is organized as follows. We illustrate in 
Sec. [n] our theoretical framework and in Sec. IIIII the 
model describing the boards-depletant mixture. Sec. IV 
is devoted to the results, whereas in Sec. [Vjwe draw our 
conclusions. 
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FIG. 1. The 6 independent orientations of a boardlike par- 
ticle within the restricted orientations (Zwanzig) model. 



to the so-called Zwanzig model: the only allowed ori- 
entations are those with the main particle axes aligned 
along the axes of a fixed reference frame [54] . Within this 
model a boardlike particle can take the 6 orientations de- 
picted in Fig. [T] and the orientation distribution function 
(ODF) is a 6-dimensional vector xp with components tpi 
(i — 1, ...,6), subject to the normalization condition 



(1) 



Being interested in the low-density phase behavior of 
the system, where the stable phases are expected to be 
homogeneous in space, we can neglect spatial modula- 
tions in the single-particle density. Under these condi- 
tions, the intrinsic free energy reads [T7] 



N 



ln(77) + y^-0»ln(V'i) 



N 



(2) 



where /3 = {kgT)^^, is the Boltzmann constant and 
rj — Nv/V is the packing fraction. A closed expression 
for the excess free energy Tcxc in terms of the ODF is not 
known in general, but for short-range pairwise additive 
potentials it is possible to write it as a virial series in rj. 
Let uw (r) be the interaction potential between a pair of 
particles with orientations i and i', respectively, and a 
separation r between their centers of mass. At second 
virial order the excess free energy reads 



II. SECOND- VIRIAL DENSITY FUNCTIONAL 
THEORY WITH RESTRICTED ORIENTATIONS 

We consider a system of N boardlike particles with di- 
mensions Ixwxt {I > w > t and particle volume v = Iwt) 
in a box of volume V at temperature T. Accounting for 
the orientational degrees of freedom at the single-particle 
level requires a numerically demanding description based 
on 3 Euler angles. In order to circumvent this problem 
while keeping the essential physics of the system, we turn 



TV = 2^ E (3) 

where the second- virial coefficients Eui /2 are given by 

Ew =- f drfu'ir), (4) 
Jv 

with the Mayer function 
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fu'{r) = exp[-/3u,,'(r)] 



(5) 



More refined approximations than Eq. ([s]) have been re- 
cently developed in order to take into account higher- 
order virial terms, which would give a quantitatively 
more reliable description [T^]. On the other hand, trun- 
cating the virial series at second order allows for an ap- 
preciable simplification of the mathematics and the nu- 
merics involved, while retaining the essential physics. 

Let us indicate with Xi the main axis (I, w or t) of 
a particle with orientation i along the x axis of a fixed 
reference frame, and similarly with Yi and Zi. Within 
the Zwanzig model each of the 6 independent orienta- 
tions of a particle can be identified by {Xi^ Yi, Zi), which 
is one of the 6 permutations of the three elements I, w 
and t. With these definitions one can write the interac- 
tion potential between two identical boardlike particles, 
modeled as hard cuboids, as 
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if |x| < iX,+X,,)/2 
and \y\ < {Y,+Y,,)/2 
and \z\ < {Z, + Z,,)/2; 

otherwise, 



(6) 



from which explicit expressions for En' in terms of I, w 
and t directly follow through Eqs. (|4])-([5]). 

At second-virial order the equilibrium ODF at fixed 
packing fraction 77 is obtained by minimizing the free en- 
ergy of Eqs. ([2|-([3]) with respect to ■0, subject to the 
normalization condition Eq. Q. In practice, the min- 
imization problem is performed by numerically solving 
the system of 6 non-linear Euler-Lagrange equations 



■01 — Cexp 



EEii'ipi, 



(7) 



with the proportionality constant C determined by Eq. 
([1]). The symmetry of the solution of Eq. ([7| allows to 
identify the stable homogeneous phase. When 0^ = 1/6 
for every i = 1, 6, the phase is isotropic (/), whereas in 
the opposite case, when all the ipi assume different values, 
the ODF describes a biaxial nematic (Nb) phase. When 
the system is characterized by the presence of a single axis 
of symmetry (uniaxial nematic phase), the coefficients 
0i are coupled two-by-two. Let us suppose this axis of 
symmetry to be the vertical axis of Fig. [T] In this case, 
we distinguish between prolate uniaxial nematic phase 
(A^+), when the most likely configurations of Fig. Il|are 
(1) and (4), and oblate uniaxial nematic phase (N^), 
when the most likely configurations are (3) and (6). 

In the treatment described so far, we assume the sys- 
tem to be homogeneous in space. In order to estimate 
the limit of validity of this assumption, we adopt bifur- 
cation theory to calculate the minimum packing fraction 



fj, beyond which homogeneous phases are unstable with 
respect to smectic states |58j. The mathematical details 
regarding the application of bifurcation theory to the free 
energy of Eqs. can be found in the Supplemental 

Material of Ref. 17 , and we report here only the final 
result. Let us indicate with Q^i^igx) the function 



QiiH^lx) = - VV'i0i' / dr fi.i,{r)exp{~iqxx), (8) 



where is the ODF of the equilibrium homogeneous 
phase at packing fraction ry, and analogously for Q^i) (qy) 

and Qii'iqz)- The bifurcation packing fraction fj^ for 
smectic fluctuations along the x axis is found as the 
minimum packing fraction at which the 6x6 matrix 

with entries Qii) (qx) has an eigenvalue 1 for some qx- 
Therefore, the smectic bifurcation packing fraction is 
fj = min[fjx,fjy,fiz)- As a final remark, it is impor- 
tant to notice that the present bifurcation analysis al- 
lows only to predict when homogeneous phases are unsta- 
ble with respect to one-dimensional modulations in the 
single-particle density. Therefore, nothing ensures the 
corresponding stable inhomogeneous phase to be charac- 
terized by one-dimensional (smectic), rather than two- 
(columnar) or three-dimensional (crystal) positional or- 
dering. 



III. EFFECTIVE DEPLETION INTERACTION 

Our aim is to study the infiuence of a depletant on the 
phase behavior of a system of boardlike particles. Hence, 
the system described in Sec. |TT] is modified by the ad- 
dition of a second species of particles (the depletant), 
modeled as cubes with dimensions d x d x d. The binary 
mixture of boardlike particles and depletant is assumed 
to be in equilibrium with a reservoir of depletant particles 
at fixed fugacity zd = &xp{l3fJ,D)/^%, where /i^ is the 
chemical potential of the depletant and Ad its thermal 
wavelength. Following the pioneering approaches to the 
topic 32 , 33^, we neglect interactions between depletants, 
in which case the fugacity zd coincides with the density 
in the reservoir. The ideal-depletant assumption is 
justified a posteriori by the low packing fractions ud d^ 
considered. Modeling the depletant with cubic particles 
appears to be rather unrealistic, especially if compared to 
typical polymeric depletants, usually treated as spheres. 
However, we claim that our choice contains the essential 
features of the physical phenomenon, while considerably 
simplifying the mathematics that follows. In the next 
section we show that the peculiar phase behavior of our 
system is due to the asphericity of the depletion volume, 
which, in turns, is a consequence of the asphericity of 
boardlike particles. Therefore, we do not expect the spe- 
cific shape of the depletion region (cuboidal for cubic de- 
pletant, spherocuboidal for spherical depletant) to play 
a major role in our results. Moreover, the relative differ- 
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ence between cuboidal and spherocuboidal depletion vol- 
ume for the values of the particles dimensions considered 
here amounts to few percentage points. The interactions 
in the mixture are given by the cuboid-cuboid potential 
Eq. ^ between boardlike particles, and by the cuboid- 
cube potential between boardlike particles and depletant, 
given by 



/3w,(r) = < 



oo 







if |x| < {X, + d)/2 
and |y| < (y. +d)/2 
and \z\ < {Z, + d)/2; 

otherwise, 



(9) 



which explicitly depends on the orientation i of the 
boardlike particle. 

At fixed fugacity Z]j the configurational entropy of the 
depletant is maximized when the total depletion volume, 
i.e. the region of space forbidden to the depletant due 
to the presence of boardlike particles, is minimized. As 
a consequence, an effective attraction between boardlike 
particles appears. Such a depletion interaction can be 
explicitly calculated by integrating out the depletant de- 
grees of freedom and must be expressed in general as a 
sum of many-body interaction terms [551 156] . For the 
sake of simplicity, we include the effect of the depletant 
by considering only the effective two-body interaction po- 
tential, while neglecting higher order terms. The effective 
pairwise depletion potential wu' (r) between cuboids with 
orientations i and i', respectively, and center-to-center 
separation r is given by [57] 



Pwu'ir) = -nDVu'ir), (10) 
with Vii' (r) the overlap volume of the depletion regions. 



if\x\> {2d + X^+Xi,)/2 

and \y\ > {2d + Y, + Y,,)/2 
and \z\ > {2d+ Zi + Zi,)/2; 

otherwise. 



Here aI^"* is defined as 



\ (^) \{y) \ (^) 
\i' \i' \i' 



(11) 



d+^^-|a;| if |x| > 

and |x| < (d + ^^-^^i' 

d + min{X,,X,,) if \x\ < 



(12) 

and analogous definitions hold for X^J^) (y) and A-^,''(z). 

Let us indicate with a tilde the properties obtained by 
adding the effective two-body depletion potential ww (r) 
to the cuboid-cuboid potential Wj,;/(r). The Mayer func- 
tion Eq. ([5| becomes 



fii'{r) = exp[-/3uii/(r) + noVwir)] - 1. 



(13) 



The phase behavior of this effective one-component sys- 
tem can then be calculated by following the prescriptions 
of Sec. [n] with the function fw (r) substituted by fwir). 
Unfortunately, the expression of Vw (r) given in Eq. ( 11 ) 



does not allow for an analytical calculation of the inte- 
grals Eii' and Q^^/^ (Qx) in Eqs. Q and However, an 
analytical expression can be obtained by truncating the 
Taylor series of the Mayer function Eq. (13) in njjVii'lr), 



h'{r) = fu'ir) + ^ -^{Vu'{r)y''e^p[-(3uu'{r)]. 



(14) 

By inserting Eq. (14) into Eq. ([4|, one obtains for the 
effective excluded-volume coefficients 



-'—'1.1,' -'—'11.' 



dviy^i' (r)) " exp [-Pu^ (r)] , 



(15) 

where the integrals of the r.h.s. can now be solved an- 
alytically for every m. Similar considerations hold for 



We verified by com- 



the functions Q'^^^) {qx) of Eq 
parison with exact numerical calculations of the effective 
excluded-volume coefficients that quantitative agreement 
can be obtained by truncating the series of Eq. (15 1 at 
fifth order in for all n^, considered in this paper. For 
consistency, the Taylor expansion in of the functions 
(5-f/'(9a;) is truncated at the same order. 



IV. RESULTS 

The framework developed in Sec. |III| allows to deter- 
mine the effective excluded-volume coefficients Eai of a 
system of cuboidal I x w x t particles due to the pres- 
ence of a cubic d x d x d depletant at fugacity zd (and 
density ud — zd). The phase behavior of this effective 
one-component system of boardlike particles is then an- 
alyzed by applying the theory described in Sec. |TTj 

It is readily understood from Eq. ( 15 ) that adding 
the depletion attraction Eq. ( 10 ) to the cuboid-cuboid 
pairwise potential Uii'{r) gives rise to a monotonic de- 
crease of the coefficients Eu' with njj. This effect is 



depicted in Fig. [2j where we report the 6 independent 
values of the matrix elements Eu, corresponding to the 6 
two-particle configurations (1,1), (1,2), (1,3), (1,4), (1,5) 
and (1,6) (cf. Fig. [T]), as a function of the reservoir 
depletant concentration njj. In order to allow for a com- 
parison with previous experimental [2S] and theoretical 
jlTj work on the subject, the aspect ratios are chosen as 
l/t = 9.3 and w/t = 3.0, while for the cubic depletant 
we set d/t = 1.0. At njj = the 6 excluded-volume ma- 
trix elements are positive definite, but with increasing 
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FIG. 2. Effective excluded-volume coefficients Eu (in units 
of boardlike particle volume v — Iwt) for the 6 independent 
orientational configurations of a pair of boardlike particles 
in the Zwanzig model as a function of the depletant num- 
ber density no- If ere the boardlike particles have dimensions 
l/t = 9.3, w/t = 3.0 and are in contact with a reservoir of 
ideal cubic depletants with side d/t — 1.0 and number den- 
sity no- 



FIG. 3. Ratio between the second-virial coefficients corre- 
sponding to the two-particle configurations (1,4) and (1,5) of 
Fig. [T]as a function of the depletant density no for boardlike 
particles with w/t = 3.0 and l/t = 8.7 {v = l/w — w/t = 
-0.1), l/t = 9.0 {u = 0.0), l/t = 9.3 {u = 0.1), and l/t = 10.0 
{v = 0.33). The solid circles highlight the value of the deple- 
tant density n^, defined by the condition E14, — Ei^. 



their value decreases until becoming negative (see Eu in 
Fig. [1]). Such a behavior is well known from the study 
of systems of spherically symmetric particles with short- 
range attractive potentials, where one can define a tem- 
perature at which the second virial coefficient changes its 
sign ("Boyle temperature"). The change in sign of the 
second-virial coefficient is related to a tendency of the 
system to develop a gas-liquid phase separation. Also in 
the present case, where the role of the (inverse) temper- 
ature is played by the depletant density nu, this change 
in sign can indicate a tendency towards a phase sepa- 
ration between two homogeneous phases. On the other 
hand, when the dimension of the depletant is sufficiently 
small, one expects the gas-liquid phase separation to be 
metastable with respect to a broad gas-solid coexistence 
[13 H?! [55J [S^ . As we ignore the stability of inhomo- 
geneous phases like smectic, columnar or crystal states, 
in the present work we limit our investigations to values 
of n£) small enough as to guarantee a positive value of 
all the effective excluded-volume matrix elements, and to 
avoid strong tendency towards a broad phase separation. 

Although the monotonic decrease with no is a fea- 
ture of all the 6 effective excluded-volume coefficients 
Eli, their rate of change is not the same. Let us fo- 
cus on the coefficients corresponding to the two-particle 
configurations (1,4) and (1,5). In absence of depletant 
{no = 0), Ei5 = Ei5 is shghtly bigger than Eu = Eu, 
but its first derivative at njj > is smaller. As a con- 
sequence, there exists a value of the depletant density 
njj, such that £^14 < Ei^ for njj < rt^, and Eu > Ei^ 
for no > n*jj (see inset of Fig. [2|. This fact, which we 
will show to have deep consequences for the phase be- 
havior of the system, is more clearly represented by the 



red solid curve of Fig. [sj representing the ratio Eu/Ex<i 
as a function oi no. In the same plot, we report the ra- 
tio Eu/Ei^i relative to boards-cubes mixtures with fixed 
ratio w/t ~ 3.0 and d/t = 1.0, but different values of 
l/t {l/t — 8.7, 9.0 and 10.0). In the four cases we ob- 
serve a monotonically increasing dependence of Eu/Eiq 
on no, which implies that the existence of njj, defined 

by the condition £'14 = £15, is determined by the value 
of Eu/ Ei^, which in turn depends only on I, w and t. In 
other words, n*jj exists only if Eu < £-15, with 71^ = 
if £14 — £15. On the contrary, if £14 > £15 one has 
£14 > £15 independently of the depletant density njy. 

Before addressing the physical consequences of the ex- 
istence of the density n^, it is worth seeing how the rel- 
ative value of the excluded volume coefficients £14 and 
£15 determines the phase behavior of boardlike particles 
in absence of depletant (i.e. njj =0). It is well-known 
that monodisperse hard boardlike particles are expected 
to undergo an IN transition, where the particles spon- 
taneously break the orientational symmetry by aligning 
along common directions in space [5 . The nematic phase 
emerging from the / can be i) uniaxial prolate N-^. with 
alignment of the long axis I; ii) uniaxial oblate with 
alignment of the short axis t; iii) biaxial Nb with align- 
ment of the three axes of the particle. Following Onsager 
[T], the origin of this phase transition can be understood 
by considering that orientational ordering determines an 
increase in excluded-volume entropy, which compensates 
the decrease in orientational entropy. Therefore at fixed 
orientational entropy, when £14 < £15 the phase will 
be thermodynamically favored over the , the opposite 
being the case when £14 > £15 (cf. Fig. [I]). In the 
intermediate situation, when £^14 = £15, the system un- 
dergoes instead a direct second-order INb transition. By 
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explicitly calculating E14 and E15 in terms of I, w and t, 
and defining a shape parameter v = l/w — w/t, one can 
show that 




< 1 <^ ly > 0, 
= 1 ^1^ = 0, 
> 1 4^ ly <0. 



(16) 



This is consistent with Straley's result that a system of 
boardlike particles undergoes i) a first-order IN^ transi- 
tion if > 0; ii) a first-order IN^ transition if < 0; iii) 
a second-order INb transition if = [B]. 

The relation between the excluded-volume coefficients 
and the character of the IN transition can be gener- 
alized to the case of boardlike particles immersed in a 
depletant, provided that the coefficients Ew are substi- 
tuted with Eiii . According to this interpretation, one is 
lead to identify the density n^, defined by the condition 
Ei4 — i?i5, with a Landau critical point at which the 
system undergoes a direct second-order INb transition. 
Since the Nb stability is generally due to a balanced com- 
petition between rodlike and platelike behavior, the crit- 
ical depletant density n*j-, is expected to divide the phase 
diagram into two distinct regions, one where the stable 
uniaxial nematic phase is prolate (A^+, corresponding to 
rodlike behavior), the others where the stable uniaxial 
nematic phase is oblate {N-, corresponding to platelike 
behavior). 

This picture is confirmed by the {ri,n]j) phase dia- 
grams of Fig. |4] describing the phase behavior of board- 
like particles with dimensions w/t — 3.0 and (a) l/t = 8.7 
[v = -0.1), (b) l/t = 9.0 {v = 0.0), (c) l/t = 9.3 
(1/ = 0.1), and (d) l/t = 10.0 {ly = 0.33) immersed in 
a cubic depletant with sides d/t — 1.0 and number den- 
sity njj- As a general feature, at packing 77 w 0.2 — 0.3 
the system undergoes a phase transition from an / phase 
(green region) to N+ (red regions), N- (blue regions) 
or Nb (yellow regions) states. The first-order charac- 
ter of the JiV_|_ and /iV_ transitions is not visible on the 
scale of Fig. |4j The dotted lines indicate the limit of 
stability of the homogeneous phases with respect to one- 
dimensional (smectic, Sm) fluctuations along the long 
axis / (red dotted lines) or along the short axis t (blue 
dotted lines). Therefore, the smectic bifurcation analy- 
sis confirms that inhomogeneous phases (white regions) 
preempt nematic states at sufficiently high packing frac- 
tions. In particular, the higher the depletant density n/j, 
the lower the stability of homogeneous phases with re- 
spect to inhomogeneous one. This result is in agreement 
with previous studies on the phase behavior of hard rods 
interacting via an attractive depletion potential [45l . 
In the latter case, in fact, the coexistence regions between 
different phases increase with the depletant density lead- 
ing, at sufficiently high njj, to a wide isotropic-crystal 
coexistence and a consequent disappearance of the liquid 
crystal phases. We expect similar phenomena at deple- 
tant concentrations higher than those considered here. 



but a description beyond the second virial order would 
be needed in this case. 

As deduced by the analysis of Fig. [3j in the case of 
"platelike" boards with j/ < (Fig. |4]ja)) the absence of 
a Landau critical point implies that the / phase under- 
goes a transition towards a iV_ phase for every value of 
nz). If we instead consider a system of boardlike particles 
with iy = (Fig. gb)), it is well known that in absence 
of depletant a direct second-order INb transition is ex- 
pected. In our picture, this corresponds to the presence of 
a critical depletant density at n*jj = 0, beyond which N_ 
states appear at intermediate packing between the / and 
the Nb phase. More interesting conclusions can be drawn 
when considering "rodlike" boards with > 0, in which 
case the critical depletant density n*jy assumes non-zero 
values (Figs. |4]jc) and (d)). If the pure system at no = 
is expected to develop an IN-^- transition, the attraction 
induced by the depletant reduces the N+ stability until 
determining a direct INb transition &i njj = n^. Sur- 
prisingly, at even higher depletant densities {no > n*£i), 
the stable uniaxial nematic phase in between / and Nb 
has oblate ordering N_, in sharp contrast with the be- 
havior of the pure rodlike boards system. Moreover, the 
phase diagrams of Figs. |4]^c) and (d) suggest that, when 
dealing with boardlike particles with z/ > 0, setting the 
depletant density at values close to allows to select 
regions of the phase diagram with relatively high Nb 
stability. This is possible also when the regime of Nb 
stability of the pure boardlike particle system is small 
(Fig. Qc)) or even absent (Fig. ^d)). 

A relevant feature of the present analysis is that a crit- 
ical depletant density n*jj exists only for slightly elon- 
gated, or rodlike, boards > 0). On the contrary, no 
Landau critical point is predicted when v < 0, in which 
case for every value oi no the system develops a first- 
order IN- transition, typical of platelike particles. This 
fact can be interpreted in the following terms. At low 
enough depletant density {ub < njj) the role of the de- 
pletant is weak and the isotropic-nematic phase transi- 
tion is driven by the gain in boardlike particles' excluded- 
volume entropy, leading to N^ (N^) ordering if 1/ > 
{ly < 0). On the other hand, at high enough depletant 
density (jib > n^) the thermodynamically more favored 
states are those maximizing the depletant entropy, i.e. 
states where the overall depletion volume is minimized. 
It appears clear that, at fixed boardlike particles' orien- 
tational entropy and independently of the sign of v, iV_ 
rather than N^ ordering tends to maximize the overlap 
between the depletion regions of single boards. There- 
fore, when ly > the Landau critical point at n*jy appears 
as a result of a competition between the excluded-volume 
entropy of boardlike particles and depletant. Instead, 
when v < this competition does not happen since both 
entropies are maximized by N^ states, and thus no crit- 
ical depletant density exists. 

One could wonder how the phase behavior of our sys- 
tem changes by varying its relevant parameters, that is, 
the dimensions of the boardlike particles Z, w, t and the 
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FIG. 4. Phase diagrams of boardlike particles with aspect ratios w/t = 3.0 and (a) l/t = 8.7 {y = —0.1), (b) l/t — 9.0 
= 0.0), (c) l/t = 9.3 {v = 0.1), (d) l/t = 10.0 {v = 0.33) in contact with a reservoir of cubic depletant with side length 
d/t — 1.0 aX number density nn- The diagrams feature isotropic (/, green regions), prolate {N+, red regions) and oblate (A''_, 
oblate regions) uniaxial and biaxial {Nb, yellow regions) nematic phases. The black circles highlight the Landau critical points, 
whereas the dotted lines indicate the limit of stability of nematic phases with respect to smectic (Sm) fluctuations along the 
long (red dotted line) and short (blue dotted line) particle axis respectively. 



size of the depletant d. The phase diagrams of Fig. |4] 
highlight the monotonic increasing dependence of n*jy on 
I/, when the dimensions of the depletant d/t and the as- 
pect ratios of one of the two boards (here w/t) are fixed. 
By numerically solving the equation En = i?i5, we in- 
vestigate also the role of the depletant dimension on the 
critical depletant density. In Fig. [5]we report the critical 
depletant density n*j^ as a function of the depletant side 
d/t for three boardlike particle dimensions with differ- 
ent aspect ratios l/t and w/t, but same shape parameter 
!/ = l/w — w/t = 0.1. As a general trend, by increasing 
the size of the depletant, the critical depletant number 
density n*jj decreases, in accordance with the intuitive no- 
tion that the smaller the depletant, the more one needs 
to establish enough depletion attraction. Moreover, at 
fixed u the critical depletion density decreases most for 
the more extreme aspect ratios of the particles. In other 
words, the bigger the aspect ratios l/t and w/t at fixed 
V, the smaller the amount of depletant needed to reach 



n\). If instead of the number density one considers the 
critical depletant packing fraction n^d'^, one sees that 
this quantity is an increasing function of d (inset of Fig. 
[5]). Therefore, one expects the ideal depletant approxi- 



mation (cf. Sec. Ill) to be increasingly reliable in the 



limit of small depletant. 

In practical situations, one could be interested in esti- 
mating the critical depletant density njj , which is defined 
as the solution of the non-linear equation _Ei4 = Ei^ and 
which has then to be calculated numerically. If njj is 
sufficiently small, one can obtain an approximate expres- 
sion for this quantity by linearizing both sides of equation 
i?i4 = i?i5 in The approximate critical depletant 

density is given by the following expression 

. ^ 2[t{l + wf~l{w + tf]{l-t)-^ 

2{lt-w^)d^ + w{lw + tw-2lt)d^' ^ ' 

and it can be compared (dotted lines of Fig. [s]) with the 



8 



0.03 



0.1 



0.01 



0.001 : 




FIG. 5. Critical depletant density njj as a function of the 
side d of the cubic depletant for different boardlike particles 
with same shape parameter v = Q.l: I /t — 9.3 and w/t = 
3.0 (solid red line); l/t = 16.4 and w/t = 4.0 (solid green 
line); l/t = 25.5 and w/t = 5.0 (solid blue line). The dashed 
lines represent the approximate analytical dependence given 
by Eq. (171. The inset illustrates the same data in terms of 



the critical depletant packing fraction n^d 



numerical calculation (solid lines), showing good overall 
agreement, which improves the larger the depletant side 
d/t. 



V. CONCLUSIONS 

In the present paper we investigate for the first time 
the effect of a short-range depletion-induced attraction 
on the liquid-crystal phase behavior of boardlike parti- 
cles. To this aim, we make use of classical density func- 
tional theory truncated at second virial order, and adopt 
the Zwanzig model for the description of the orientational 
degrees of freedom. In close analogy with the Asakura- 
Oosawa-Vrij model for mixtures of spheres, by neglecting 
interactions between the cubic depletant particles we can 
explicitely calculate the effective two-body attractive de- 
pletion potential between boardlike particles. 

We predict that in systems of slightly elongated board- 
like particles (i' > 0), there exists a critical depletant den- 
sity at which the uniaxial nematic phase is substituted by 
a direct second-order transition from the isotropic to a bi- 
axial nematic phase. At higher depletant concentrations. 



a large region of oblate uniaxial nematic ordering devel- 
ops, rendering the system of attractive rodlike boards be- 
having like a system of hardly repulsive platelike boards. 
The origin of this phenomenon is due to two competing 
mechanisms: the maximization of the boardlike particle 
entropy, favoring 7V+ ordering, and the maximization of 
the depletant entropy, favoring N_ ordering. 

The phase behavior described in this work shares 
many similarities with our findings in Ref. [T7j, where 
we showed that size-polydispersity in a system of hard 
boardlike particles with the same shape and different 
volume induces the appearance of a Landau tetracrit- 
ical point at a specific system composition. This fact 
is related to a competition between prolate and oblate 
ordering, which in turn is realized only when the board- 
like particles are slightly elongated. In the light of our 
present findings, we suggest that this N+-N^ competi- 
tion and the corresponding emergence of a Landau tetra- 
critical point can be understood in terms of a depletion 
effect. More specifically, when size-polydispersity be- 
comes relevant, N- rather than iV+ ordering determines 
the higher total entropy due to the minimization of the 
overall depletion regions of the big particles with respect 
to the smaller ones. In further analogy with the boards- 
depletant mixture, no such competition is predicted for 
platelike boards and consequently no tetracritical point 
appears in this case. 

Besides furnishing an explanation for the results of Ref. 
[17] , we suggest that manipulating the attraction induced 
by a depletant, e.g. a non-adsorbing polymer, furnishes 
an original and effective way to control the phase behav- 
ior of boardlike particles, allowing to stabilize prolate and 
oblate uniaxial and biaxial nematic states. Moreover, the 
depletant density is expected to be an easy experimental 
control parameter. 
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